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1. Introduction 

Students tend to regard the elongated "S-shaped" [H El [8] logistic curve of pop- 
ulation dynamics (fig. I) as somewhat exotic. It is typically derived by applying 
the method of partial fractions to a separable differential equation [I] 1^1 131 ISl HI 13 
m m [TOl [11] . My purpose here is to show how the logistic curve may be derived 
more directly as a simple consequence of the more familiar differential equation 
model for exponential decay, and that the curve itself is nothing more than a famil- 
iar friend in disguise. The disguise is removed by abandoning our fixation on the 
reference point (0, Pq), representing the initial population at time zero, in favour of 
a much more natural choice. This illustrates an important principle, namely that 
one should always adapt the coordinates to the problem at hand. In this case, a 
great deal is simplified by relocating the origin more appropriately. 



2. Background 

Textbooks (see eg. fT][ll[21[Sl[71[Hl[ni[ini[n]) typically begin the discussion of 
population growth with the exponential model 

1 dP 

p^ = k =^ P(i) = Poe'=*, t>0, (2.1) 

in which the relative growth rate fc is a positive constant, representing say the av- 
erage birth rate. Since unbounded growth is unrealistic, more sophisticated models 
take into account limited resources for reproduction. The logistic model, proposed 
by the Belgian mathematical biologist Pierre F. Verhulst in 1838 [T], replaces the 
constant relative growth rate in (j2.ip with a relative growth rate that decreases 
linearly as a function of P: 

1 dP ( P\ 

--^fc(l--j, .>0, 0<Po<Af. (2.2) 

The constant M represents the maximum sustainable population beyond which P 
cannot increase. The dimensionless factor 1 — P/M in (|2.2p serves to diminish the 
relative growth rate from k down to zero as the population increases from its initial 
level Po to M. 

Although one can solve (|2.2p as a Bernoulli differential equation by making 
the substitution P ~ 1/y [3], for the most part texts treat (|2.2p as a separable 
differential equation to be solved by the method of partial fractions [U [21 [31 [51 
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H [HI [9l [TOl [H] . Either way, one obtains, after some algebraic simplifications, the 
solution 

nt) = ^-^—^, where R,^^i-Il. (2.3) 
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Figure 1. Logistic Curve with M = 10, Po = l, k = 0.7 



3. Approach Via Exponential Decay 



Suppose instead of counting individuals, we count niches, viewing Af as the 
maximum number of niches the ecosystem can support, and P as the number of 
niches currently occupied. Dimensional analysis suggests that instead of studying 
P, we consider 

M-P , , 

i?=— (3.1) 

the dimensionless ratio of available or vacant niches to niches currently occupied. 
From ()2.2|) . it readily follows that 



'R'dt 



R{t) = Roe 



-kt 



(3.2) 
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where Rq is as in (12. 3p . Expressing p.2|) in terms of P, one arrives again at the 
sohition ([O]) . 

For the instructor who would hke to discuss logistic growth but would prefer 
to bypass separable differential equations and partial fractions, it may be desir- 
able to "cut to the chase" by introducing the logistic model via (|3.ip and (|3.2[) 
rather than via (|2.2|) . In other words, although the fact that R satisfies the dif- 
ferential equation (|3.2p follows from its definition (|3.ip and the logistic differential 
equation (j2.2p . one could equally well dispense with V2. 2\) and make LS. S\} an as- 
sumption of the model, thereby proceeding more quickly to the solution. Of course, 
there are pedagogical advantages to either approach. One aspect the approach 
via ()3.ip and (13. 2p we are proposing has in its favour is that logistic growth can be 
introduced in the standard section on exponential growth and decay, with no loss 
in continuity and without any additional background. 

That R decreases at a rate proportional to itself, i.e. satisfies the differential 
equation (|3.2p . is intuitively plausible. Initially we think of P being much smaller 
than M, so that R is much larger than 1 and many niches are available relative 
to the number currently occupied (a high niche vacancy rate). We should expect 
any species to take advantage of such a hospitable climate for reproduction, and 
hence initially, R should decrease rapidly as P increases. However, as the number 
of vacancies decreases, {P gets close to M, R gets close to zero) there are relatively 
few available niches remaining. In such an inhospitable climate, we should expect 
reproduction and hence further growth to be difficult, and accordingly, R should 
decrease much more slowly. These considerations should be sufficient to motivate 
the introduction of logistic growth via (13. ip and ()3.2p to any calculus or differential 
equations class. 

4. Removing the Disguise 

From the viewpoint of an individual of the species attempting to reproduce, one 
should expect a qualitative change in the hospitality of the ecosystem near R = 1, 
given the considerations of the previous paragraph. Motivated by these consider- 
ations, we refer to an ecosystem as being hospitable or inhospitable according to 
whether R is greater or less than 1. From (|3.2p . the transition from hospitable to 
inhospitable occurs when 

^ - T-o := ^ logi?o ^ ^ log { ^'^pf° ^ . ^=1, P^^M. (4.1) 

It is well-known that this is precisely the time at which P is increasing most rapidly, 
as can be seen by completing the square in (|2.2p : 

dP _ Mk k / mV 
'dt ~ ~ ^ M \ YJ ' 

Because of the distinguished nature of the point (ro,M/2) it seems more sensible 
to measure time from tq than from zero. Certainly t = is completely arbitrary 
from the viewpoint of the species, having more to do with whatever external forces 
(desire, opportunity, availability of funding etc.) conspired to allow the biologist 
or census taker to obtain an initial field count than any essential features of the 
system. Therefore, we consider 

Q(r) F(to + t) - iikf. 
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Pit) - ^^-H , ^0 > 0. 



where r measures time from tq and hence may be positive or negative. From (12. 3|) . 
we have 

n( \ - ^ ^ 

Since i?oe~'^^° = 1, this simpUfies to 

Q[t) = ^'^ , -E = E( 1^^^ = iMtanh(i/tT) , (4.2) 
or in other words, 

P{t) = iM (l + tanh (ifc(t - To))) , (4.3) 

where tq is given by (|4.ip . Thus, the mysterious "S-shaped" [21 [HI [5] logistic curve is 
nothing more than a translate of our old and familiar friend, the hyperbolic tangent. 

5. Addendum 

If < Po < , then oo > i?o > 0. The boundary cases Pq = and Pq = M 
correspond to Rq — oo, P{t) = and Rq — 0, P{t) = M, respectively. To complete 
the analysis of logistic growth, it is necessary to consider what happens when Pq 
lies outside the closed interval [0, Af], i.e. Rq < 0. The solution (|2.3p is valid for 
such Rq, but (j4.3|) was predicated on the assumption i?o > in the definition of tq. 
Putting So — —Rq, we have from (|2.3p that 

M 

1 - Sqc 
In this case, we define 

To := T log'S'o so that e'^^" = S'o - ^° (5.1) 
k Pq 

A calculation analogous to (|4.2p reveals that 

Hir) P(ro + r) - iAf = iMcoth (i/cr) , 

or 

P{t) = \M (1 + COth {\k{t - Tq))) , (5.2) 

where now tq is given by (|5.1|) . 

If Po > M, then -1 < Po < 0, < 5o = -Po < 1, and -oo <tq <0. Therefore, 
for t > we are on the upper arch of the hyperbolic cotangent, with population 
decreasing exponentially to M as i ^ cxo (fig. 2). In the less biologically meaningful 
case Pq < 0, we have — 1 > Pq > — oo, 1 < 5*0 = — Pq < oo and < tq < oo. As t 
increases from zero to tq, the rightmost portion of the lower arch of the hyperbolic 
cotangent is traversed, sending the population to minus infinity. The asymptote is 
then crossed and we skip over to the upper arch, the population reverting to its 
behaviour in the previous case. 
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Figure 2. Logistic Curve with M = 10, Pq = 30, fc = 0.7 



